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CAN A SATELLITE GALAXY MERGER EXPLAIN THE ACTIVE PAST OF THE GALACTIC CENTER? 
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ABSTRACT 

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating 
that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion- 
powered activity - this paints a rather different picture from the GC as we understand it today. We examine a 
possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky- 
Way which began at the redshift of z sa 10 and ended few million years ago with a merger of the Galactic 
supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite. 

Subject headings: galaxies: interactions — Galaxy: center — Galaxy: kinematics and dynamics — Galaxy: 
nucleus 



L INTRODUCTION 

There is mounting observational evidence that the epoch 
that ended several million years ago was marked by an un- 
usual level of activity in the Galactic Center (GC). This is 
remarkable given that at the current epoch, the GC is best 
characterized by the quiescent and underluminous nature of 
Sgr A* . The picture of the GC as once powerful nucleus has 
begun to emerge from circumstantial observational evidence, 
most recently strengthened by a discovery of the "Fermi bub- 
bles", a pair of giant gamma-ray emitting bubbles that extend 
nearly 10 kpc north and south of the GC (Dobler et al.|20I0 



Su et al. 2010). The Fermi Bubbles were almost certainly 
triggered by an abrupt episode, and current explanations in- 
clude a past accretion event onto the superm assive black hole 
(SMBH; |Su et al.|20I0[ [Zubovas et al.|20I I|l, AGN jet s ( |Guol 



|& Mathews|20Il t, a nuclear starburst (! Su et al.|20I0|l, and a 
sequence of star capture events in the last ~ 10 Myr ( Cheng et 
al.|20II I. The period of increased gamma-ray activity is con- 
sistent with the finding that until several hundred years ago 
Sgr A* was orders of magnitude more X-ray luminous than it 
is today, as indicated by the echo in the fluorescent Fe K line 
emission detected in t he direction of the molecular clouds in 
the vicinity of S gr A* CIn ui et al |2009| JPonti et al |20I0[[l^ 
rier et al.|20I0i. Although we cannot be certain that the cur- 



rent quiescence of the GC is unusual, it appears clear that the 
GC experienced an active phase as recently as a few hundred 
years ago. 

The GC is also a hotbed of star formation containing the 
three most massive young star clusters in the Galaxy: the Cen- 
tral cluster, the Arches cluster, and the Quintuplet cluster (see 
[Figer 2008 , for a review). The three clusters are similar in 
many respects. Each cluster contains ^ IQ*Mq in stars and 
has central stellar mass density that exceeds those measured 
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in most globular clusters. While our current understanding 
of massive star and star cluster formation is incomplete, it is 
plausible that these clusters are all characterized by the star 
formation event within the past 2-7 Myr that resulted in for- 
mation of more massive st ars (above 100 Mq) than anywher e 
else in the Galaxy (K rabbe et al.|I995[|Pauniard et al.|2006| ). 
It is possible that the three clusters have a common origin and 
that they have formed as a consequence of a single event that 
triggered the flow of the copious amo unts of gas into th e cen- 
tral ^ 50 pc in the Galaxy (though see Stolte et al.[ 2008 for a 
scenario in which the Arches cluster forms at the intersection 
of X1/X2 gas orbits in the inner Galaxy) 

On even smaller scales, the existence of massive and young 
stars in the Central cluster, well within the central parsec, is 
especially puzzling. Although a view that the young stel- 
lar population formed in situ seems to be prevailing, it is 
still a mystery how the starforming clouds that gave birth 
to this population escap ed tidal disruption by the central 
SMBH UPaumard et~aL] [2006). A clue in favor of the in- 
situ formation is that most young, massive stars seem to in- 
habit one or more disk-like structures, pointing to their birth 
in a dense accretion disk ^Bartko et al.||2010| l. Star forma- 
tion in a gaseous disk also provides a natural explanation for 
the cuspy distribution of the young stars. In order for star- 
forming clumps to withstand tidal forces in the inner parsec 
of the GC their densities need to be in excess of 10"cm~^, 
at least five orders of magnitude higher than the average den- 
sity of molecular clouds in the GC (Figer et al._20 00). Such 
densiti es can only be achieved through highly compressive 
events ( Figer|20()8 1 and it is plausible that both the inflow of 
large amounts of gas into the GC and its shocking and com- 
pression have been caused by a common culprit. Both phe- 
nomena are found to arise as consequences of galactic merg- 
ers (Noguchi 1988; Barnes & Hernquist 1991, 1992 1996| 
iMihos& Hemquist,1996, .Hopkins & Quataei't,2010| ) making 
this a possibility worth examining. 

While the young stellar distribution appe ars to be 
cuspy (Genzel et al. 2003; Paumard et al. 2006; Buchholz et| 
al. 2009; Do et al. 2009; B artko et al.|20I0| l, there seems to be 
a distinct lack of old stars. This evidence is based on num- 
ber counts of spectroscopically identified old stars brighter 
than magnitude K = 15.5 within the sub-parsec region about 
Sgr A* . The best fits of the inner density profile for old popu- 
lation seem to favor power-laws with slopes of 7 < 1 and even 



allow the possibility of a core with 7 < 0, with the stellar den 
sity decreasing toward the cent er (Buchholz et al.|[2009 Do 



et al.||2009l [Bartko et aL]|2010 |). At this stage, the evidence 



for the deficit of old stars is compelling, however there are 
still significant uncertainties in their density profile: the pop- 
ulation of stars based on which this inference has been made 
are luminous late-type giants that comprise only a small frac- 
tion of the underlying stellar density of the old population. 
Regardless of the precise slope, however, the distribution of 
the old population is contrasted by the steeply rising density 
distribution of young stars. 
Possible mechanisms that could create a core in the distri- 



bution of old stars have been discussed by Merritt ( 2010^ and 
include destruction of stars on orbits that pass close to the 
SMBH in a triaxial nucleus, inhibited star formation near the 
SMBH at the time when the old population was formed, and 
ejection of stars by a massive black hole binary. In light of the 
other evidence that points to a discrete event in the recent his- 
tory of the GC, it is interesting to revisit the latter mechanism. 
In giant elliptical galaxies, the existence of cores is of- 
ten attrib uted to ejection of stars by an inspiralling b inary 
SMBH UFaber et al.|1997[|Milosavljevic & Me rritf2001 ) and 
possibly due to gravitational wave recoil after binary coa- 
lescence ([Boylan-Kolchin et al. 2004, Gualandris & Merntt] 
2008]!. The prediction of these models is that the central steT 



lar mass deficit (traced by the stellar light) is proportional to 
the mass of the central black hole, Mjef oc M,. This correla- 
tion is interesting in view of the observed dichotomy between 
ellipticals with cores and those with the extra central light: 
core light deficit was found to correlate closely with M, and 
stellar velocity dispersion a, in agreement with the theoreti- 
cal predictions (" Milosavljevic & Merritt|200T Gualandris & 
Merritt 2008), however, the extra light does not (Kormendy 
S Bender 2009 ). An expl anation of these phenomena offerea 
by |Kormendy & Bender| (200 9 ) is that the extra light ellip- 
ticals were made in wet mergers with starbursts, where stars 
formed from gas leftover after the merger, whereas core el- 
lipticals were created in dry mergers. In galaxies with excess 
light, the newly formed population of stars fills the core left in 
the distribution of the older population to form a steep cusp, 
thus giving rise to the different distributions of old and young 
stars, not unlike that indicated in the GC. 

While studies of the light excess and deficit in elliptical 
galaxies focus on major mergers, this seems less relevant for a 
disk-dominated system like the Mil ky Way, which may h ave 
never experienced a major merger ( [Gilmore et al.|[2002] l. A 
minor merger of the SMBH with an intermediate mass black 
hole (IMBH) however cannot be ruled out. The presence of 
an IMBH in the Galactic center has been previously consid- 
ered asa^ossibleveliicl£fOTdeliyery of young stars into the 
GC ( Hansen & Milosavljevic'2003'), a mechanism for creation 
of hypervel ocity stars (Baumgardt et al. 2006 ), and growth of 
the SMBH (Por tegies Zwart et al.|2006[ ). Indeed, a possibiHty 
that an IMBH with mass < IO'^Mq is still lurking in the inner 
parsec of the GC cannot be totally excluded based on obs er- 
vations (H ansen & Milosavljevic|2003[IGenzel et al.|2010 i. 

In light of the new observational evidence, which supports 
the notion that few to 10 Myr ago was a special period in 
the life of the Sgr A*, as indicated by the relatively recent 
episode of star formation and increased energy output, we re- 
visit the possibility that a minor merger could have triggered 
this epoch of enhanced activity. We suggest that the cumula- 
tive observational evidence favors the minor merger hypoth- 



esis relative to the scenarios that propose a steady state evo- 
lution or passive relaxation of the GC region. A theoretical 
scenario for such minor merger is discussed in § 2 and discus- 
sion is presented in § 3. 

2. MILKY WAY - SATELLITE MERGER SCENARIO 

Here we examine the viability of the following scenario: 
at redshift ^10, a primordial satellite galaxy with a central 
IMBH begins to merge with a young Milky Way. As the 
satellite sinks toward the GC under the influence of dynam- 
ical friction it is tidally stripped and its orbit gradually de- 
cays toward_flie_^4il^ Way disk plane (Quinn & Goodman] 
[1986) [Callegari et~aL][20Tl] l. The satelhte perturbs previ- 
ously stabl e gas clouds in the inner Milky Way disk, driving 
gas inflow ( Noguchi" 19881 [Barnes & Hernq uist|[T99T] [T9961 

'"'1 



Hopkins & Quataert 2010} and compressing TEe gas to den- 



sities exceeding those ne cessary for massive star formation 
near t he GC ( ,Mihos& Hernquist||1996t [Hopkins & Quataert 



2010|l. The satellite galaxy is expected to be largely disrupted 



by the time it reaches the GC, leaving the IMBH spiraling 
in a dense gaseous and stellar environment. In the context 
of this scenario we hypothesize that the IMBH reached the 
central parsec on the order of ^10 Myr ago. A fraction of 
perturbed gas t hat did not form stars acc retes onto the Milky 
Way's SMBH (Hopkins & Quataert'2010il, injecting massive 
amounts of energy into t he sur rounding medium and giving 
rise to the Fermi bubbles (! Suet al.|20 10; Zubovas e tal.|2011[ ). 
Once gravitationally bound, the IMBH-SMBH binary orbit 
tightens via three-body interactions with surrounding stellar 
backgrou nd, scouring th e old stellar population to form a cen- 
tral core ( Merritt|2010 i. Finally, the binar y coalesces follow- 
ing the emission of gravitational radiation ([Peters & Mathewsl 

\im\ . 

In the context of this hypothetical scenario we use the new 
GC observations to constrain the initial masses of the satellite 
and Milky Way galaxies (Msat and Mmw), the mass deficit in 
the old stellar population (Mdef), the IMBH mass (Mimbh), as 
well as the amount of gas inflow into the GC triggered by the 
inspiral of the satellite galaxy. 

2.1. Properties of Merging Galaxies 

Beginning the merger at high redshift is advantageous in 
three respects. First, at this early epoch, it is reasonable to 
assume that the proto-MiUcy Way was surrounded by primor- 
dial satellite galaxies capable of h ousing a central seed black 
hole (e.g.,[Ri cotti & Gnedin 2 005ljGnedin & Rravtsov|2006 
Wise & Abel|2008MMicic et al.|2(]l II ). Second, at this stage 
in its growth, the Milky Way would have been smaller, less 
massive, and more gas-rich than it is today, thus decreasing 
the time required for the satellite to sink to the galactic center 
via dynamical friction. Finally, the orbits of infalling satel- 
lites are more radial at high redshi ft, which further shortens 
the merger timescale ( Wetzel|20TT I. 

To determine the properties of the Milky Way at this epoch, 
we assume that it grows according to the exponential halo 
model from [McBride et al.[ ( [2009| l: 



M(z) = Mo(l+z/exp(-ln2 



z/ 



(1) 



where Mq is the current halo mass, and Zf is the formation 
redshift, defined as the redshift at which the halo has grown 
to half its current mass. Adopting the properties for the 
Milky Way at z/ = 1 as Mo = 2 x lO'-M© and (3 = 0.25, the 
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Milky Way's mass at z = 10 can be estimated to be Mmw = 

In order to deliver the IMBH to the GC a mere 2-7 Myr ago, 
the proposed merger redshift of ^10 implies that the satellite 
orbit decayed over a timescale of about 1 3 Gyr At such a high 
redshift, the IMBH and satellite had very little time to evolve 
before being accreted by the Milky Way, making the pair a 



could drive ^ IQt'Mq of gas toward the center. Given that 
even a few % perturbation can destabihze the gas and drive 



"fossil" of dark ages before reio nization ( Ricotti & Gnedin 
2005 Gnedin & Kravtsov 2006[l. To constrain the merger 



history of the satellite, however, several more ingredients are 
necessary in addition to the merger time scale: the mass ra- 
tio Msat/Miviwi as well as the energy and angular momentum 
of the satellite orbit. We rely on cosmological N-body simu- 
lations to constrain the initial conditions of the orbit. These 
inform us that at the present epoch, satellites are preferentially 
accreted on very eccentric orbits, with a distribution peak at 



about e = .85 (Benson 2005; Wang et al. 2005 ; Zentn er et al. 
2005j[Khochfar & Burkert,2006, .Ghigna et al.,1998,,Tormen 



et al. 1997 1. At higher redshifts the satelhte orbits are char- 
acterized by even higher eccentricities, albeit, in both cases 
the distribution peaks are broad. Seemingly independent of 
redshift, a typical satellite is accreted with a total velocity 
Vtot = 1.15Vvir (Vvir IS the circular velocity at the Viri al radius 
of the primary galax y), that marks it as barely bound ( |Benson| 
2005 Wetzel|20lT i. Motivated by these results in the context 
of our model we select an orbit that has Viot = 1.1 5 Kir at the 
Virial radius of the primary and an eccentricity of 0.97, con- 
sistent with expectations for the eccentricity distribution peak 
at z = 10 (Wetzel 201 1 ). Given these assumptions, we set out 
to determine the mass of the satellite galaxy. 

We calculate the orbital decay for a range of satellite masses 
placed at the Virial radius with an energy and angular mo- 
mentum defined above. We adopt a static model of the 
Milky Way that incl udes a central SMBH, Miyamoto-Nagai 
thin disk (Miya moto & Nagai||197 5), a sph erical Hernquist 
bulge (Hernquist 1 1990| l, and an NFW halo ( jNavaro, Frenk,| 
[& White|1997| ). To mimic a young Milky Way, we scale the 
mass and size of the Milky Way halo from z = 0toz=10 
(with the NFW halo concentration c = 10) in such way that 
the ratio of the mass of baryons to the halo is preserved. 
We integrate the orbits inc luding dynamical fric tion as de- 
scribed by Chandrasekhar ( [Chandrasekhar [r943 1 and adjust 
the Coulomb logarithm as mass is lost from the satellite, 
A = 1 + (Msat/Mpri^ary)^- We include tidal stripping to the 
Roche radius of the satellite and the effect of mass loss as the 
satellite passes through the disk plane. However, we neglect 
the stellar component of the satellite, since the baryon content 
of such low mas s satellites is relatively uncertain, but likely 
to be very s mall ( |Gnedin|2000l [Simon & Geha|2007| |Ricotti| 
|et al.|2008| l. We find that the rnost likely culprit is a satellite 
with a mass of Msat ~ 5 x lOP Mq. Modeling the satellite dark 
matter profile as an NFW halo, its corre sponding concentra- 
tion parameter at this redshift is about 50 ([Bullock et al.|2001 1 



Our estimate indicates that by the time the satellite has 
reached the inner 100 pc, it will have lost ^97% of its mass, 
with only ^ 2 x 1O^M0 remaining. A ^ IQ^Mq satellite is 
massive enough to disturb the gas in the Milky Way that is 
accumulated in a ring at the Inner Lindblad Resonance (ILR) 
at the radius ^150 pc ( Stark et al. 2004|l. Using disk sta- 
bility theory developed by 'Papaloizou et al. (2007), we es- 
timate that the orbiting satellite remnant of this mass could 
excite a ^ 50% perturbation on the gas ring at the ILR, which 



inflow ("Barnes & Hernquist"! 99 1| [1996[ [Mihos& Hernquist 
1996, Hopkins & Quataert 2010]i7the stripiped satelhte/IMBH 
pair would be massive enough to deplete the ILR of gas as it 
plunges past it. 

2.2. Late-Type Stellar Mass Deficit and IMBH Mass 

If the core in the distribution of old stars at the GC ( [Buch-| 
'holz et al. 2009[[Do et al.[2009| l was scoured out by an IMBH- 
SMBH binary, the amount of stellar mass missing from the 
GC can be used to constrain the mass ratio of the black hole 
binary dMilosavljevic & Meifitt1[2001[[Gualandris & Merritt| 
2008| l. To determine this mass deficit we compare the stellar 
distribution inferred from observations with that expected for 
a dynamically relaxed system without a core. In terms of the 
number density of the late-type stellar population the core can 
be represented by a broken power law 



«/(r) = «o 



'"o 



1 + 



'"o 



(7,-7)/q 



(2) 



withno = 0.21 pc"-', tq = 0.21 pc, 7= 1.8, 7,- = -1.0, and a = 4 
(iMerritt 2010). We adopt this description in our analysis but 
note that in presence of strong mass segregation the slopes 



can be steeper ( Alexander & Hopman 2009 ; Preto & Amaro 
'Seoane[ [2010| [Amaro-Seoane & Preto||2011| ). The observed 
distribution of stars outside of the 0.21 pc co re radius is con- 
sisten t with the Bahcall-Wolf profile (ex r"' ^ Bah call& Wolf 
[1976| l of a relaxed system as it would have existed prior to 
scattering by the IMBH-SMBH binary. We model the initial 
stellar cusp by extending the Bahcall-Wolf profile to smaller 
radii: 

ni{r) = nJy\ . (3) 

Assuming that the mass density profiles before and after the 
creation of the core are given by equations ([3]l and ([2| re- 
spectively, we calculate the mass deficit as the integrated dif- 
ference between the initial and final (observed) profiles. We 
choose the mass density normalization of 1 .0 ±0.5 x 10^ M0 
inside the inner parsec, as determined by Schodel et al. 
(2009), and obtain Mdef « 2 x lO^M©. 

N-body merger simulations studying the relationship be- 
tween the ratio of total stellar mass ejected to binary mass, 
Mief/{M\ +M2), and binary mass ratio, q = M1/M2, have not 
yet been carried out for the mass deficit calculated here. In 
what follows we use a semi-analytic formalism describing the 
interaction of massive bla ck hole binaries with their stellar en- 
vironment ( Sesana et al.[ [2008ll to place the upper and lower 
limits on the mass of the IMBH based on Mdef inferred from 
observations. 

It has been sho wn by numerical s imulations (Baumgardt et| 
laL 200 6j,Matsubayashi et al.|20 07t and semianalytic models 
( Sesana et al.l 



that an IMBH inspiralling in a stellar 
cusp surrounding a central SMBH starts to efficiently eject 
stars at a separation aq, where the stellar mass enclosed in the 
IMBH orbit is of the order of 2M2. The ejection of bound 
stars causes an IMBH orbital decay of a factor of « 10, ex- 
cavating a core of radius kq w 2ao in the central stellar cusp, 



resulting in a mass deficit of the order of 3M2 (see Sesana et 
al. 2008 for details). Such orbital decay is in general insuf 
ficient to bring the IMBH in the efficient gravitational wave 
(GW) emission regime, unless its eccentricity grows to > 0.9 



during the shrinking process. It is also the case in this picture 
that the mass of the inspiralling IMBH inferred for a given 
mass deficit strongly depends on the eccentricity evolution of 
its orbit. In what follows, we consider both, the high and low 
orbital eccentricity scenario and use them to place a bound on 
the plausible range of IMBH masses. 

If the eccentricity grows efficiently, the IMBH depletes the 
central cusp, forms a core of a size « 2ao, and mer ges due to 
GW emission on a time scale of only 1-10 Myr (Sesana et 



long lived (> 1 Gyr) binary configuration at separation < rp 



|al.|2008 ). For a stellar distribution described by an isothermal 
sphere outside of the radius of influence of the SMBH, oo = 
Iq^/^pc. Adopting the core radius of kq = Iuq = 0.2 Ipc, we 
find q = 0.02, and an upper limit on the mass of the IMBH, 
M2 = 8 X 1O"^M0. In this case, the mass evacuat ed from the 
stellar cusp by the IMBH is of the order of 3M2 ( Sesana et al.' 
[2OO8J , i.e., « 2.5 X 1O^M0, consistent with the stellar mass 
deficit measurement in the GC. 

Alternatively, if the IMBH eccentricity does not grow sig- 
nificantly during the bound cusp erosion, further scattering of 
stars replenishing the binary loss cone is needed in order to 
evolve from separation of ao to the GW regime. Therefore, 
a circular orbit regime can be used to establish a lower limit 
on the mass of the IMBH, for a given mass deficit indicated 
by observations. In this case we assume that in this case both 
ro and Mjef created in the cusp erosion phase are small (we 
justify this assumption below). In this scenario, the final ro 
and Mdef are reached as a consequence of the diffusion of the 
stars from the edge of the small core into the loss cone of the 
binary. The ejection s of each star carry away an energy of the 
order (3/2)G/i/fl (Q uinlan||1996| , where fi = M1M2/M. We 
compute Mdef by imposmg: 



3 Gfj, 



dM^ 



GM1M2 1 
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to get 



Ml +M2 fly 

3 fl/ 



(4) 



(5) 



where fli is the hardening radius of the binary (radius at which 
the scattering of unbound stars becomes effective) and flf is 
the separation at which the GW emission bec omes efficient. 
Using equations (19) and (20) in Sesana ( 2010| to express Oi 
and flf, it follows that. 



M1+M2, 5000^/5 

A/* = In — -, 

3 F(e)i/5 ' 



(6) 



where F(e) = (l-e2)-7/2(i+ 73/24^2 _,_ 37/96 g4)^ Assuming 
for the purpose of this estimate that the binary remains circu- 
lar throughout its evolution and imposing Mdef = 2 x lO^M©, 
we find q = 5 X 10"^ and a lower limit on the mass of the 
IMBH M2 = 2 X IO^Mq n An IMBH of such mass, would 
excavate a core of w O.Olpc, causing a mass deficit of ~ 
6 X IO-'M© in the bound scattering phase and thus, justifying 
our earlier assumptions. Note that in the circular orbit sce- 
nario the time scale for the inspiral of the IMBH towards the 
GW regime is determined by the unknown rate of diffusion of 
the stars into the loss cone of the binary. Hence, depending 
on the time scale of relaxation processes this process could 
in principle lead to the IMBH-SMBH binary "hangup", i.e., a 

' Both numerical simulations and semi-analytic models however suggest 
that the eccentricity in the cusp erosion phase grows to > 0.9, in which case 
F(e)> 1000and?>5x lO"^ i.e.,M2 >2x 1O''M0. 



f 



tantamount to the classical "final parsec" problem (Begeknan 
|etal.|1980| l. 

It is however unlikely that the binary will stall in our spe- 
cific case. The galactic center in this phase will be described 
by a strongly perturbed, non-axisymmetric potential which al 
lows stars to scatter into the loss cone efficiently (Merritt 8i\ 
[Poonl|2004l IBerczik et all[2006l |Perets & Alexander||2008 
Moreover, the orbit will occur in a relatively gas-rich environ- 
ment, which can further aid the decay of the binary ( [Escala et| 
ial. 2005, Dotti et al. 2007 ; Cuadra et al.; 200 9 ). Finally, any 
extra stars brought in by the satellite would help the binary 
decay (see ! Miller,2002| . Even under the assumption of a cir- 
cular orbit, an efficient coalescence on a timescale of 10 Myr 
is therefore plausible. 

This analysis suggests that the observed mass deficit and 
core size are consistent with the IMBH mass in the range 
2 X IO^^Mq < M2 < 8 X IO^Mq, whereas the eccentricity evo- 
lution found in N-body simulations and semi-analytic models 
favor M2 > 1O"^M0. Within this range, the time scale for the 
IMBH to create a core and merge with the SMBH can be as 
short as few Myr On the other hand, a possibility that an 
IMBH may be still be lurking in the GC is also allowed. We 
discuss the consequences of the latter scenario in the context 
of the observational constraints on the presence of a second 
black hole in the Galactic center in §[3] 

It is useful to consider whether a satellite galaxy with an 
initial mass of Msat -- 5 x lO^M© can host a > IO'^Mq IMBH. 
While there are no observational constraints for galaxies or 
black holes of this mass range, there are three leading theo- 
ries for IMBH formation at high redshift: 'direct collapse' of 
metal-free, low angular momentum gas into a 10^- 10^ Mq 
black hole (Loeb & Rasio„1994, Begelman et al.|l2008]l , an 
unstable su permassive star that collapses into a 10"^ - 10^ Mq 
black hole (|Colgate"1967t [Quinlan & Shapiro [[19871 [Baum-| 
o 19991, or a Population III star, which 
lO^Mp 



garte & Shapirc 



leave behind seed black holes of 10- 



would 
between red- 



shift 30-12 (Madau & Rees.2dOT][Bromm et al.|2002) [Wise] 
[& Abel[[2008t [Clark et al.[[2'oi 1| ). Even if the IMBH in our 
satellite started as a low mass seed formed from a Pop III star, 
it is plausible that it would reach an IMBH mass proposed 
here. This is because the satellite galaxies at z = 10 can have 
a significant fraction of their mas s is gas, with a baryon f rac- 
tion comparable to the universal ( Ricotti & Gnediii|2005 1. In 
such satellite galaxy, it would require less than one percent 
of the gas to accrete onto a low mass seed to form the IMBH 

Note that the massive seeds produced in a direct collapse 
typically favor more massive halos than the one we have pro- 
posed as our culprit. This is because metal-free gas collapses 
most efficiently in halos with T^j^ > 10"* K, corres ponding to 
Mvii- > lO^M0[(l+z)/lO]^/2 (Bromm & Loeb'2003 ). In the 
context of the merger hypothesis choosing a slightly more 
massive satellite would push the accretion redshift closer to 
the present day, and as long as the resulting satellite merger is 
still a minor one, this does not significantly affect the outcome 
of our scenario. 

2.3. Inflow of Gas and Gamma-ray Bubbles 



As noted in § 2. 1 the inspiral of a satelhte galaxy can cause 
the inflow of a si gnificant amount of gas towards the cen- 
ter of the G alaxy (Noguchi'1988'; 'Barne s & Hemquist[T996 



[Mihos& Hernquist 1996, Cox et al. 2008| l. One fraction of 
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this gas could have given rise to the star formation in the 
Central, Arches, and Quintuplet clusters, which marked the 
epoch between 2-7 Myr ago in the central 50 pc of the Milky 
Way. All three clusters contain some of the most massiv e stars 
in the Galaxy and have inferred masses of ^ 10'^ Mq (Figer 
|2008|l . Assuming a "standard " star formation efficiency of 



10% ( jRownd & Young||1999| l, it follows that the amount of 
gas necessary to produce the stellar population of the three 
clusters is few x 10^ Mq. Note that a sequence of strongly 
compressional events during the satellite-Milky Way merger 
coul d have given rise to a higher efficiency of star forma- 
tion ( |Di Matteo et al.|2007| l, in which case the estimated mass 
of the gas represents an upper limit. 

In this merger scenario, the remainder of the perturbed gas 
that did not form stars would be channeled towards the cen- 
tral parsec, and the fraction that is accreted into the SMBH 
could drive the energetic outburst of several Myr ago. The 
far-IR and millimeter observations indicate that ^ 10"* M© of 
the molecular gas continues to reside in the circu mnuclear 
disk within the central ~' 1.5 pc of the Galaxy (see Genzel 
|et al.|[2010] for review and references therein). The maxi- 
mum amount of the remnant molecular gas that has not been 
accreted onto the SMBH can also be estimated based on its 
expected gravitational effect on the orbits of the stars resid- 
ing within the inner 0.5 pc. In this case, the requirement for 
stability of the stellar disc over its lifetime of 6 Myr poses 
a constraint on th e mass of the molecular torus of < 10^ Mq 
dSubr et al.|2009l). 



On the other hand, the recent discovery of the two large 
gamma-ray bubbles extending from the GC above and be- 
low the galactic plane are compelling evidence of a relatively 
recent period of intense activity in the now quiet GC. The 
gamma-ray bubbles exhibit several striking properties: they 
are perpendicular and symmetric with respect to the plane 
of the Galaxy, have nearly uniform gamma-ray brightness 
across the bubbles, and well defined sharp edges ( [Dobler e"t| 
laLI2010l|Su et al. 2010). The gamma-ray emission from the 
bubbles is characterized by the hard energy spectrum and is 
most likely to originate from the inverse Compton scattering 
of the interstellar radiation field on the cosmic ray electrons 
- the same population of electrons deemed responsible for 
the diffuse synchrotron microwave radiation detected by the 
WMAP (Fi nkbeiner 2004, Dobler & Finkbe iner 2008). The 
sharp edges of the Fermi bubbles are a lso trac ed by th e X-ray 
arcs discovered in the ROSAT maps ( [Snowden et al.||1997| l, 
suggested to be the remnants of shock fronts created by the 
expanding bubbles ( |Su et al.|20l"0l|Guo & Mathews|201 1] . 

The morphology, energetics, and emission properties or the 
Fermi bubbles favor the explanation that bubbles were created 
in a strong episode of energy injection in the GC in the last 
^10 Myr that followed an accr etion event onto the SM BH 
( |Su et al.|20"T0| . Simulations by |Guo & Mathews] ( [20TT] ) in- 
dicate that the bubbles could have been formed by a pair of 
bipolar jets that released a total energy of 1 - 8 x 10^^ erg over 
the course of ^0.1-0.5 Myr between 1 and 2 Myr ago. This 
explanation for the Fermi bubbles implies that ^ W^Mq of 
material must have been accreted onto the SMBH at nearly 
the Eddington rate, assuming the accretion efficie ncy of 10% 
(Shakura & Sunyaev 1973 Davis & Laor 2011"). Based on 
the range of models explored by |(juo & Mathews] (J20TT]) it 
is possible to estimate that the amount of mass processed in 
such jets (i.e., the mass of the gas that fills the jet cavities) is 
as small as 30Mq and as large as 3 x 10^ Mq. 



This estimate, together with the gas that formed stars, the 
gas accreted onto the SMBH and the gas processed by the 
jets allows us to put a constraint on the total gas inflow into 
the central ^ 50 pc of the Galaxy of < IQ^ Mq, consistent 
with the amount expected from the perturbation analysis of 
the stability of the ILR gas in the Milky Way. 

3. DISCUSSION 

3.1. How rare are satellite merger events? 

We propose that the timeline began about 13 Gyr ago, when 
the proto Milky Way accreted a small satellite dark matter 
halo at the time when their halos were physically smaller and 
closer in mass. The satellite orbit decayed slowly and only 
reached the GC a few million years ago, after having been 
stripped of most its mass. The thinness of the Milky Way disk 
has ofte n been used as an argu ment against a re cent minor 
merger (fQuinn, Hernquist, & Fullagar 1993; Sellwood, Nel- 



son, & Tremaine 1998; Velazquez & White 1999); however, 



the proposed satellite is so minor, particularly by the time the 
orbit decays to 10 kpc, th at the thi n disk cou ld have survived 
unscathed (Toth & Ostriker' 1992', "Walk er et al.|1996[ [Taylor] 
& Babul 2001 ; Hopkins et al. 2008 2009). Using cosmologi- 
cal N-body simulations to derive merger rates for M ilky Way 
mass halos as a function o f mass ratio and redshift ( |Sinha &| 
|Holley-BockeImann|201 l) l, we can estimate how often these 
minor accretion events occur Over the course of the lifetime 
of a typical Milky Way halo, it has suffered 10-30 satellite 
accretions in the Ms^t = 10^- lO^M© range. Only between 10 
and 20% of these can sink to the galactic center by the present 
day, separated by a time scale of ^ few Gyr. 

Figure 3. 1 illustrates a realization of the current distribution 
of accrete d satellites in the mass range Msat = 10^ - 1O^M0. 
We used |Sinha & Ho lley-Boc kelmann| (2011) to define the 
number of accreted halos at a given redshift for this mass 
range and set the energy and angular momentum distribu- 
tions at each redshift using the expressions 7-9 in (Wetzel] 
201 l|l. As in section 2, we integrated orbits of the satellites 



from the accretion epoch to the present day, scaling the M ilky 
Way mass and size to the redshift of accretion. Figure |3.1| 
shows the inner 40 kpc of the current-day Milky Way; ap- 
proximately 2/3 of the accreted satellites are at separations 
larger than 40 kpc, while 15% have already reached the GC 
and merged with the SMBH. We estimate the surface bright- 
ness of the satellites assuming that the baryons are confined to 
a radius within the dark matter halo ten times smaller than the 
satellite virial radius. We infer the initial star fraction from Rj3 
|cotti & Gnedin (2005 ) and assume that the mass-to-light ratio 
is 3 for the bound stars that remain after the tidal stripping. 

While it is clear that not all of the small satellites can rich 
the GC, what fraction does is a question of some subtlety. 
Galaxy merger timescales cited in the literature, particularly 
for the small mass ratios considered here, span a wide range. 
The key to the uncertainties is the treatment of dynamical fric- 
tion: most semi-analytic works, including this one, rely on the 
dynamical friction formalism as described by Chandrasekhar 
(1943) but change the Coulomb loga rithm to account for in- 
homogeneous or anisotropic systems ( Penarrubia et al.|2004[ 
Just & Penarrubia 2005 ), or to include mass loss ( [Taylor &| 
' Babuli |2001 ; Velazquez & White[[1999] ). This approach has 



been shown to un derestimate the decay time in pure dark 
matter simulations (Colpi et al.[[1999] Boylan-Kolchin et al. 



2008 1. On the other hand, the presence of gas can dramat- 



ically decrease the orbital decay time of a satellite (Ostriker] 
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Fig. 1 . — Distribution of accreted low mass satellites at the present day. The inner 40 kpc region of the Milky Way disk is shown in greyscale with cun'ent 
accreted satellite positions overplotted. The color maps to the surface brightness, and the relative size corresponds to the tidal radius of the satellite. Note that the 
circle size is not to scale and that none of these satellites would be observable above the background. 



1999 Sanchez-Salcedo & Brandenburg|1999 1 - thus, making 
the Chandrasekhar formula a significant overestimate. Com- 
pounding the issue, linear perturbation theory and limited N- 
body experiments indicate that resonant heating caused by or- 
bits in the satellite galaxy that are commensurate with the or- 
bit of the satellite about the GC can enhance mass loss and 
can change the angular momentum of the orbit in non-trivial 
ways (Weinberg 1997; Choi et al. 2009). Although the Milky 
Way has likely accreted dozens of these small satellites, it 
is uncertain how often they reached the galactic center It 
is however plausible that the GC has experienced several of 
these accretion events spread over its lifespan. Naturally this 
uncertainly in the dynamical friction timescale translates into 
an uncertainly in the mass of our culprit satellite, as well. Fu- 
ture ultra-high resolution hydrodynamic nbody simulations of 
satellite accretion in a growing Milky Way will help to pin 
down minor merger timescales, which will help to further 
constrain the range of plausible sulprit satellite masses and 
accretion redshifts. 

3.2. Hypervelocity stars and stellar core 

While the properties of the newly formed stars and per- 
turbed gas were dictated by the accreted satellite, the IMBH 
was responsible for carving out the old stellar population. As 
a gravitationally bound binary IMBH-SMBH formed and de- 
cayed, it scoured out 2 x lO^M© of the relaxed old and ini- 



tially cuspy stellar population. Many of these stars could 
have been ejected from the GC as hypervelocity stars (HVSs; 
[Baumgardt e t al. 2006), though most may simply have re- 
ceived enough energy to traverse the inner parsec. Simula- 
tions of IMBH-SMBH binaries in stellar environments indi- 
cate that HVSs are created in a short burst which lasts only a 
few Myr in case of a ~ lO'^M© IMBH ('Baumgardt etal.'2006| 
Sesana et al.|2008 i. In the context of this picture we predict 
that this event created ~ 10^ hyper velocity stars that, if they 



were ejected at about 1000 kms"' ( Baumgardt et al.||2006 l 
ought to lie a few kpc from the GC today. It is worth not- 
ing that about a dozen of HVSs observed in the Galactic halo 
thus far have travel times that span 60-240 Myr and appear to 
be consistent with a continuous ejection model ([EFo wn 2008] 
Brown et al.|[2009l ITilhch et al.][2009| [Irrgang et al.,.201 



and not with the IM BH-SMBH bmary picture ( |Brown|2008| 
|Sesana et al.|2008 V Along similar lines, the spatial and veloc- 
ity distribution of the observed HVSs seem to be inconsistent 
with a IMBH-SMBH sHngshot origin (Sesana et al. 2007). It 
is however possible that the observed HVSs have been ejected 
from the GC before the IMBH reached the central stellar cusp 
in the Galaxy and are thus not associated with the binary. 

The large size of the observed GC core, ro = 0.21 pc, could 
be seen as a challenge to any scenario involving 3-body scat- 
tering, since state of the art high resolution direct N-body sim- 
ulations that modeled the ejection of hypervelocity stars from 



Active Past of the Milky Way 



a SMBH-IMBH binary in the gala ctic center never ge nerated 
a core larger than 0.02 parsecs (Baumgardt et al. 2006| l. How- 
ever, there are several effects that could conspire to cause the 
simulated core size t o be a lower limit. First, the mass of the 
simulated SMBH in|Baumgardt et al.|(2006) is 3 x IO'^Mq, 



For 10'* M0 IMBH the black hole merger can give rise to 
a modest recoil velocity of about 80ms"', assuming that the 
IMBH is not spinning rap idly. The recoil velo city mag nitude 



in this case scales as ex q (jCampanelli et al 



2007; Baker etal. 



which would eject fewer stars than somewhat more massive 
Milky Way SMBH. Second, the density profile was sharply 
curtailed by a factor of (1 +r'') in order to minimize the num- 
ber of stars far from the SMBH; this makes the spatial dis- 
tribution of stars in the simulated nuclear star cluster more 
centrally peaked relative to that in the GC, which can also re- 
sult in a smaller core. In general, though, it is important to 
note that the size of the scoured core is a property that sensi- 
tively depends on the density, the eccentricity, and kinematic 
structure of the GC or on assumptions in the model used to 
represent it. 

3.3. Has IMBH-SMBH binary merged? 

We now return to the question whether the IMBH-SMBH 
binary has already merged or whether t he I MBH could still 
be lurking in the GC. As discussed in § |2.2| the N-body and 
semi-analytic modeling of the GC favor the evolutionary sce- 
narios in which the inspiral and coalescence of the SMBH 
with a M2 > IQ^Mq IMBH is relatively efficient. Moreover, 
there is currently no empirical evidence for a second black 
hole in the central parsec. In order to be consistent with the 
observations, the IMBH present in the GC would have to have 
a mass ^ lO-'- lO^^M© and be eit her very close (< lO '^pc) 
or at > O.lpc from the SMBH (see |Genzel et al.||2010 l, for a 
review of observational constraints). An IMBH in this mass 
range that reaches a separation of lO'^'pc would merge with 
the SMBH in less than 10 Myr due to the emission of GWs, 
thus significantly restricting the amount of parameter space 
where the IMBH and SMBH can exist in a long lived binary 
configuration. Nevertheless, given the uncertainties in the bi- 
nary mass ratio, eccentricity, and the structure of the initial 
stellar cusp, the presence of an IMBH in the GC cannot be 
entirely ruled out at this point. 

If on the other hand, the IMBH and SMBH coalesced sev- 
eral million years ago, one possible signature of this event 
could be a SMBH recoil caused by the asymmetric emission 
of GWs (Peres 1962; Bekenstein 1973). Current astromet- 
ric observations of the reflex motion of the SMBH put strong 
constraints on the allowed recoil velocity; the SMBH cannot 
have velocity with respect to the Central cluster larger than 
3.5kms~' ( within \cr er ror), at the distance of the GC (Yelda 
et al.| |2010). Similarly, |Reid & Brunthaler] ( |2004 ) constrain 
the peculiar motion of Sgr A* in the plane of the Galaxy 
to -18 ±7kms"' and perpendicular to the Galactic plane to 
-0.4±0.9kms"', where quoted uncertainties are \a errors. 
There is a however a caveat with respect to the interpretation 
of the SMBH reflex motion: if the reference frame in which 
the reflex motion is measured is based on the nearby gas and 
stars bound to the SMBH, the resulting relative velocity of 
the SMBH will be zero because in this case, the stars and the 
gas move together with the SMBH as long as their orbital ve- 
locity is higher than the that of the re flex motion. The r adio 
and near-infrared reference frames in [ Yelda et aT] ( |2010[ ) are 
defined based on the nearby stars orbiting around the SMBH 
and are thus a s ubject to this caveat. The measurement of 
|Reid & Bruntha ler (2004) is however carried out in the refer- 
ence frame defined by the extragalactic radio sources and can 
be used to test the recoil hypothesis. 



|2008| ), thus implying that the coalescence of the SMBH with 
a slowly spinning IMBH more massive than 1.5 x 10^ M0 can 
be ruled out based on larger of the observational constraints, 
as long as damping of the recoil motion of a remnant SMBH 
is inefficient on the time scale of several million years. More 
stringent constraints on the mass of the IMBH, based on the 
motion of the SMBH perpendicular to the Galactic plane, can 
be placed given the (unknown) orientation of the orbital plane 
of the binary before the merger in addition to the binary mass 
ratio and the spin vector of the IMBH. 

3.4. Orientation of the SMBH spin axis 

The nearly perpendicular orientation of the spin axis of the 
SMBH to the Galactic disk plane, indicated by the orienta- 
tion of the observed gamma-ray bubbles and jets in simula- 
tions of Guo & Mathews (201 1 ), implies that the evolution of 
the SMBH spin has been determined by accretion from the 
Galactic gas disk rather than random accretion events with 
isotropic spatial distribution. Such events would include tidal 
disruptions of stars and giant molecular clouds triggered by 
the satellite inspiral and a merger with the satellite IMBH 
which orbital plane in principle may not be aligned with the 
plane of the Galaxy. It is thus interesting to consider whether 
a sequence of such accretion events can exhibit a cumulative 
torque on the SMBH sufficient to re-orient its spin axis, as- 
suming that before the merger with a satellite galaxy it was 
perpendicular to the Galactic plane. 

Consider first the effect of episodi c gas accretion resulting 
from multiple tidal disruption events. |Chen et aL] ( |2009||201 1 1 
show that three-body interactions between bound stars in a 
stellar cusp and a massive binary with properties similar to 
the IMBH-SMBH considered here can produce a burst of tidal 
disruptions, which for a short period of time (^ 0.1 Myr) can 
exceed the tidal disruption rate for a single massive black hole 
by two orders of magnitude, reaching N ~ 10"^ yr"' . This im- 
plies that in the process of the IMBH inspiral the SMBH could 
have disrupted ^ 10^ stars. A key element in this considera- 
tion follows from the finding by Natarajan & Pringle ( 1998} 
and Natarajan & Armitage^ (1999) that the orientation of the 
spin axis of a SMBH is very sensitive to the angular momen- 
tum of the accreted gas: namely, accretion of a mere few % in 
mass of a SMBH can exert torques that change the direction 
but not the magnitude of the spin of a black hole. Because 
each in a sequence of random accretion events imposes an in- 
finitesimal change in the orientation of the SMBH spin axis, 
collectively they can cause the spin axis to perform a random 
walk about its initial orientation. Thus, the magnitude of the 
effect scales with the number of disrupted stars and their mass 
as ^ \/Nm^. Since this is much less than few percent of Mi, 
the cumulative effect of tidal disruption events on the orienta- 
tion of the spin axis of the SMBH will be negligible. 

This conclusion is reinforced by an additional property of 
post-tidal disruption accretion disks: they are compact in size 
and usually confined to the region of a size few x r,, where 
r, « r^, (Mi/m^,y/^ is the tidal disruption radius of a star and 
r* is the stellar radius ( Rees|[T988) . Such small accretion 
disks effectively act as very short lever arms for torques act- 
ing on the spin axis of the SM BH, thus furth er reducing the 
efficiency of this process (Natarajan & Pringle, 1998 ). 
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Similar conclusion can be reached in case of the tidally dis- 
rupted molecular clouds and the gas flows that plunge towards 
the SMBH on nearly radial orbits as a consequence of the per- 
turbation by the satellite galaxy. In section § |2.3| we estimated 
that the amount of mass accreted by the SMBH is ^ 10'^ Mq. 
A modest mass, combined with the small circularization ra- 
dius of the gas accretion disk is insufficient cause a signifi- 
cant change in the SMBH spin orientation. Even "accretion" 
of a spinning IMBH is not expected to noticeably influence 
the spin orientation of the remnant SMBH. The large mass 
ratio of the binary ensures that the final contribution of the 
IMBH's spin and orbital angular momentum to the final spin 
of the SMBH is small, as long as pre-merger SMBH has a 
moderate initial spin, > few x 0.1, in terms of th e dimension- 
less spin parameter (Barausse & Rezzolla|2009) . Hence, the 
coalescence with the IMBH would not have had a significant 
effect of the SMBH spin axis orientation. 

In summary, the torques from the accretion of tidally dis- 
rupted stars, gas, and IMBH in the aftermath of the satellite 
inspiral will be insufficient to change the orientation of the 
SMBH spin axis as long as the SMBH spin is > few x 0. 1 . 
It follows that the perpendicular orientation of the spin axis 
has been set by the physical processes before the merger with 
the satellite, and most likely by the accretion of gas from the 
Galactic disk. 

4. CONCLUSIONS 

A range of theoretical arguments and observational evi- 
dence make a compelling picture for a satellite infall event 
within our GC which triggered a brief epoch of strong star- 
forming and AGN activity millions of years ago. When cou- 
pling the newest data - on the Fermi bubble and the dearth 
of old stars - to the well-established features of the GC such 
as the cuspy young star population, a timeline of the recent 
dynamical events in the galactic center emerges. 

While the case for a merger of the Milky Way with a satel- 
lite galaxy is not ironclad, it is a plausible explanation that 
naturally accounts for both the old and young stellar distribu- 
tions and the recent violent past of Sgr A*. This event may not 
be unique in the evolution of the Milky Way; indeed N-body 
simulations of the growth of Milky Way-mass galaxies sug- 
gest that the present epoch is rife with mergers of relic satellite 
galaxies wi th the galactic center, occurring at a rate of one per 
few G yr (Diem and et al.|2007[ Sinha & Holley-Bockelmann 
201 1 1. This implies that there may have been other bursts of 
hypervelocity star ejections, which can seed a population of 
intragroup stars" farther out in the halo of the Galaxy. Inter- 



estingly, we see tentative evidence in the SDSS archive for a 
potential set of very late M giants at ^ 300kpc , outside the 
virial radius of our galaxy ( [Palladino et al.|20TT i. Although a 
followup observation is needed to ensure that these intragroup 
candidates are not L dwarfs, if these do prove to be very dis- 
tant giants, they may be provide supporting evidence of a pre- 
vious minor-merger induced burst of ejected stars ^10** years 
ago. 

Along similar lines, if satellite infall induced activity is 
common, then there may be a subset of spiral galaxies which 
exhibits the signs of the recent onset of the accretion-powered 
jets. While the longer term X- and 7-ray signatures of jets 
expanding into the intragalactic medium may be too faint 
to observe in galaxies other than the Milk y Way, relatively 
bright and short lived radio-jets (^ 0. 1 Myr; |Guo & Mathews] 
2011 1 may be present in a fraction of up to ~ 10""* Milky- 



Way-Uke spirals, assuming the minor merger rate cited above. 
Some of these galaxies may be observed serendipitously, dur- 
ing the transient phase associated with the onset of a power- 
ful jet, similar to the case of the previously inactive galaxy 
J164449.3H-573451 that was recently detected by the Swift 
observ atory as a powerful source of beamed emission ( |Bur^ 
,rows et al.|20I 1 1. If it can be shown that such a sequence of 
events occurred in not so distant past in our Galaxy, it would 
forever change the paradigm of the Milky Way as an inactive 
galaxy with an underluminous central SMBH. 



We thank Cole Miller, Melvyn Davies, Jorge Cuadra and 
Rainer Schodel for insightful discussions. K.H-B., T.B., RA- 
S. and A.S. also acknowledge the hospitality of the Aspen 
Center for Physics, where the work was conceived and car- 
ried out. K.H-B. acknowledges the support of a National 
Science Foundation Career Grant AST-0847696, and a Na- 
tional Aeronautics and Space Administration Theory grant 
NNX08AG74G as well as the supercomputing sup- port of 
Vanderbilt's Advanced Center for Computation Research and 
Education, and NASA's Pleiades and Columbia clusters. Sup- 
port for TB. was provided by the National Aeronautics and 
Space Administration through Einstein Postdoctoral Fellow- 
ship Award Number PF9-00061 issued by the Chandra X- 
ray Observatory Center, which is operated by the Smith- 
sonian Astrophysical Observatory for and on behalf of the 
National Aeronautics Space Administration under contract 
NAS8-03060. Support for M.L. was provided by a National 
Science Foundation Graduate Research Fellowship. 



REFERENCES 



Alexander T., Hopman C, 2009, ApJ, 697. 1861 

Amaro-Seoane R, Preto M., 201 1, Classical and Quantum Gravity, 28, 

094017 
Bahcall, J. N. & Wolf, R. A. 1976, ApJ, 209, 214 
Baker, J. G., Hoggs, W. D., Centrella, J., Kelly, B. J., McWilliams, S. T., 

Miller, M. C, & van Meter, J. R. 2008, ApJ, 682, L29 
Barausse, E., & Rezzolla, L. 2009, ApJ, 704, L40 
Barnes, J. E. & Hemquist, L. 1991, ApJ, 370, L65 
Barnes, J. E., & Hernquist, L. 1992, ARA&A, 30, 705 
Barnes, J. E., & Hernquist, L. 1996, ApJ, 471, 115 
Bartko, H., et al. 2010, ApJ, 708, 834 
Baumgarte, T. W., & Shapiro, S. L. 1999, ApJ, 526, 941 
Baumgardt, H., Gualandris, A., & Portegies Zwart, S. 2006, MNRAS. 372, 

174 
Begelman, M. C. Blandford, R. D., & Rees, M. J. 1980, Nature, 287 
Begelman, M. C, Rossi, E. M., & Armitage, R J. 2008, MNRAS, 387, 1649 
Bekenstein, J. D. 1973, ApJ, 183, 657 



Benson. A. J. 2005, MNRAS, 358, 551 

Berczik, R, Merritt, D., Spurzem, R., & Bischof H. R 2006, ApJ, 624, 21 

Boylan-Kolchin, M., Ma, C.-R, & Quataert, E. 2004, ApJ, 613, L37 

Boylan-Kolchin, M., Ma, C.-R, & Quataert, E. 2008, MNRAS, 383, 93 

Bromm, V., & Loeb, A. 2003, ApJ, 596, 34 

Bromm, V., Coppi, R S., & Larson, R. B. 2002, ApJ, 564, 23 

Brown, W. R. 2008, arXiv:08 11.0571 

Brown, W. R., Geller, M. J., & Kenyon, S. J. 2009, ApJ, 690, 1639 

Buchholz, R. M., Schodel, R., & Eckart, A. 2009, A&A, 499, 483 

Bullock, J. S., et al. 2001, MNRAS, 321, 559 

Burrows, D. N., et al. 201 1 (arXiv:l 104.4787) 

Callegari, S., Kazantzidis, S., Mayer, L., Colpi, M., Bellovary, J. M., Quinn, 

T., & Wadsley, J. 201 1. ApJ, 729, 85 
Campanelli, M., Lousto, C, Zlochower, Y.,& Merritt, D. 2007, ApJ, 659, L5 
Chandrasekhar, S. 1943, ApJ, 97, 255 

Chen, X., Madau, P, Sesana, A., & Liu, F. K. 2009, ApJ, 697, L149 
Chen, X., Sesana, A., Madau, R, & Liu, F. K. 201 1, ApJ, 729, 13 



Active Past of the Milky Way 



9 



Cheng, K. S., Chernyshov, D. O., Dogiel, V. A., Ko, C. -., & Ip, W. -. 201 1, 

ApJ, 731,L17 
Choi, J.-H., Weinberg, M. D., & Katz, N. 2009, MNRAS, 400, 1247 
Clark, P. C, Glover, S. C. O., Klessen, R. S., & Bromm, V. 2011, ApJ, 727, 

110 
Colgate, S. A. 1967, ApJ, 150, 163 

Colpi, M., Mayer, L., & Governato, F. 1999, ApJ, 525, 720 
Cox, T. J., Jonsson, R, Somerville, R. S., Primack, J. R., & Dekel, A. 2008, 

MNRAS, 384, 386 
Cuadra, J., Armitage, P. J., Alexander, R. D., & Begelman, M. C. 2009, 

MNRAS, 393, 1423 
Davis, S. W. & Laor, A. 201 1, ApJ, 728, 98 
Diemand, J., Kuhlen, M., & Madau, P 2007, ApJ, 667, 859 
Di Matteo, P, Combes, R, Melchior, A.-L., & Semelin, B. 2007, A&A, 468, 

61 
Do, T., Ghez, A. M., Morris, M. R., Lu, J. R., Matthews, K., Yelda, S., & 

Larkin, J. 2009, ApJ, 703, 1323 
Dobler, G., Finkbeiner, D. P, Cholis, I., Slatyer, T., & Weiner, N. 2010, ApJ, 

717, 825 
Dobler, G., & Finkbeiner, D. P 2008, ApJ, 680, 1222 
Dotti, M., Colpi, M., Haardt, F, & Mayer, L. 2007, MNRAS, 379, 956 
Ernst, A., Just, A., & Spurzem, R. 2009, MNRAS, 399, 141 
Escala, A., Larson, R. B., Coppi, P S., & Mardones, D. 2005, ApJ, 630, 152 
Faber, S. M., et al. 1997, AJ, 1 14, 1771 
Figer, D. F, et al. 2000, ApJ, 533, L49 
Figer, D. F 2008 (arXiv:0803.1619) 
Finkbeiner, D. P 2004, ApJ, 614, 186 
Genzel, R., et al. 2003, ApJ, 594, 812 
Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Reviews of Modem 

Physics, 82, 3121 
Ghez, A. M., et al. 2008, ApJ, 689, 1044 
Ghigna, S., et al. 1998, MNRAS, 300, 146 
Gillessen, S., Eisenhauer, R, Trippe, S., Alexander, T., Genzel, R., Martins, 

F, & Ott, T. 2009, ApJ, 692, 1075 
Gilmore, G., Wyse, R. E G., & Norris, J. E. 2002, ApJ, 574, L39 
Gnedin, N. Y. 2000, ApJ, 542, 535 
Gendin, N. Y. & Kravtsov, A. V. 2006, ApJ, 645, 1054 
Gualandris, A., & Merritt, D. 2008, ApJ, 678, 780 
Guo, E & Mathews, W. G. 201 1, Draft 
Hansen, B. M. S., & Milosavljevic, M. 2003, ApJ, 593, L77 
Hernquist, L. 1990, ApJ, 356, 359 
Hopkins, P E et al. 2008, ApJ, 688, 757 
Hopkins, P E, Cox, T. J., Younger, J. D., & Hernquist, L. 2009, ApJ, 691, 

1168 
Hopkins, P E & Quataert, E. 2010, MNRAS, 407, 1529 
Inui, T., Koyama, K., Matsumoto, H., & Tsura, T. G. 2009, PASJ, 61, S241 
Irrgang, A., Przybilla, N., Heber, U., Nieva, M. E, & Schuh, S. 2010, ApJ, 

711, 138 
Just, A., & Peiiarrubia, J. 2005, A&A, 431, 861 
Khochfar, S., & Burkert, A. 2006, A&A, 445, 403 
Kormendy, J., & Bender, R. 2009, ApJ, 691, L142 
Krabbe, A., et al. 1995, ApJ, 447, L95 
Loeb, A., & Rasio, F A. 1994, ApJ, 432, 52 
Madau, P, & Rees, M. J. 2001, ApJ, 551, L27 
Matsubayashi, T., Makino, J., & Ebisuzaki, T. 2007, ApJ, 656, 879 
McBride, J., Fakhouri, O., & Ma, C. 2009, MNRAS, 398, 1858 
Merritt, D. 2010, ApJ, 718, 739 
Merritt, D.& Poon, M. Y. 2004, ApJ, 606, 788 
Micic, M., Holley-Bockelmann, K., & Sigurdsson, S. 2011, MNRAS, 414, 

1127 
Mihos, J. C. & Hernquist, L. 1996, ApJ, 464, 641 
Miller, M. C. 2002, ApJ, 581, 438 
Milosavljevic, M., & Merritt, D. 2001, ApJ, 563, 34 



Miyamoto, M. & Nagai, R. 1975 , PASJ, 27, 533 

Noguchi, M. 1988, A&A, 203, 259 

Natarajan, P, & Armitage, P J. 1999, MNRAS, 309, 961 

Natarajan, P, & Pringle, J. E. 1998, ApJ, 506, L97 

Navaro, J. F, Erenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Ostriker, E. C. 1999, ApJ, 513, 252 

Palladino, L., Holley-Bockelmann, K., Durrell, P., Morrison, H., Feldmeier, 

J., Wade, R. 201 1, submitted to Astronomical Journal 
Papaloizou, J. C. B., et al. 2007, Protostars and Planets V, 951, 655 
Paumard, T., et al. 2006, ApJ, 643, 1011 
Peres, A. 1962, Phys. Rev., 128, 2471 
Perets, H. B., & Alexander, T. 2008, ApJ, 677, 146 
Peiiarrubia, J., Just, A., & Kroupa, P 2004, MNRAS, 349, 747 
Peters, P C, & Mathews, J. 1963, Phys. Rev., 131, 435 
Ponti, G., Terrier, R., Goldwurm, A., Belanger, G., & Trap, G. 2010, ApJ, 

714,732 
Portegies Zwart, S. F, Baumgardt, H., McMillan, S. L. W., Makino, J., Hut, 

P, & Ebisuzaki, T. 2006, ApJ, 641, 319 
Preto M., Amaro-Seoane P, 2010, ApJ Lett., 708, L42 
Quinlan, G. D., & Shapiro, S. L. 1987, ApJ, 321, 199 
Quinlan, G. D. 1996, New Astronomy, 1, 35 
Quinn, P J. & Goodman, J. 1986, ApJ, 309, 472Q 
Quinn, P J., Hernquist, L., & Eullagar, D. P 1993, ApJ, 403, 74 
Rees, M. J. 1988, Nature, 333, 523 
Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 
Ricotti, M. & Gnedin, N. Y. 2005, ApJ, 629, 259 
Ricotti, M., Gnedin, N. Y, & Shull, J. M. 2008, ApJ, 685, 21 
Rownd, B. K. & Young, J. S. 1999, ApJ, 118, 670 
Sanchez-Salcedo, F J., & Brandenburg, A. 1999, ApJ, 522, L35 
Schodel, R., Merrit, D., & Eckart, A. 2009, A&A, 502, 91 
Sellwood, J. A., Nelson, R. W., & Tremaine, S. 1998, ApJ, 506, 590 
Sesana, A., Haardt, E, & Madau, P 2007, MNRAS, 379, 45 
Sesana, A., Haardt, E, & Madau, P 2008, ApJ, 686, 432 
Sesana, A. 2010, ApJ, 719, 851 
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Simon, J. D. & Geha, M. 2007, ApJ, 670, 313 
Sinha, M. & Holley-Bockelmann, K. 2011, (arXiv: 1103. 1675) 
Snowden, S. L., et al. 1997, ApJ, 485, 125 
Stolte, A., Ghez, A. M., Morris, M., Lu, J. R., Brandner, W., & Matthews, K. 

2008, ApJ, 675, 1278 
Su, M., Slatyer, T. R., & Finkbeiner, D. P 2010, ApJ, 724, 1044 
Subr, L. ., Schovancova, J., & Kroupa, P 2009, A&A, 496, 695 
Stark, A. A., Martin, C. L., Walsh, W. M., Xiao, K., Lane, A. P, & Walker, 

C.K. 2004, ApJ, 614, L41 
Taylor, J. E. & Babul, A. 2001, ApJ, 559, 716 
Terrier, R., et al. 2010, ApJ, 719, 143 

Tillich, A., Przybilla, N., Scholz, R.-D., & Heber, U. 2009, A&A, 507, L37 
Tormen, G., Bouchet, F R., & White, S. D. M. 1997, MNRAS, 286, 865 
Toth, G., & Ostriker, J. P 1992, ApJ, 389, 5 

van den Bosch, F C, Lewis, G. F, Lake, G., & Stadel, J. 1999, ApJ, 515, 50 
Velazquez, H. & White, S. D. M. 1999, MNRAS, 304, 254 
Walker, I. R., Mihos, J. C, & Hernquist, L. 1996, ApJ, 460, 121 
Wang, H. Y, Jing, Y P, Mao, S., & Kang, X. 2005, MNRAS, 364, 424 
Weinberg, M. D. 1997, ApJ, 478, 435 
Wetzel, A. R. 201 1, MNRAS, 412, 49 
Wise, J. H. & Abel, T. 2008, ApJ, 684, 1 
Yelda, S., Lu, J. R., Ghez, A. M., Clarkson, W., Anderson, J., Do, T, & 

Matthews, K. 2010, ApJ, 725, 331 
Zentner, A. R., Kravtsov, A. V., Gnedin, O. Y, & Klypin, A. A. 2005, ApJ, 

629,219 
Zubovas, K., King, A. R., & Nayakshin, S. 201 1, MNRAS, L274 



